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Abstract. We present a study of the formation of clustered, massive galaxies at 
large look-back times, via high resolution spectroscopic imaging of CO in the unique 
GN20 proto-cluster. The data reveal evidence for rich structure and gas dynamics in 
unprecedented detail, allowing us to image the molecular gas with a resolution of only 
1 kpc just 1.5 Gyr after the Big Bang. These state-of-the-art data give new insight 
into the detailed physical processes involved in early massive galaxy formation, and 
they provide a first glimpse of the morphological studies that will become feasible on a 
regular basis with ALMA. 



1. INTRODUCTION 



The majority of submillimeter galaxies ( SMGs: iBlain et all 120021) are believed to be 
starburst-dominated major mergers (e.g., IChapman et al. 20031) . Indeed, there is di- 
rect evidence for multiple CO component s and/or disturbed kinematics in some SMGs, 
supporting the merger pict ure (e.g. jTacconi et al.ll2008l : lEngel et al.ll2010l : llvison et al. 
201 ll : iRiechers et al.ll20lTb . However, it has recently been suggested that other mech- 
anisms that drive extreme star formation rate s may also be at pl ay. In particular, the 
cold mode accretion phenomenon ( CMA; e.g. Dekel et al. 2009a.b) has been extended 
to SM Gs by a number of authors dFardal et al.ll200ll : iFinlator et al.l 120061 : iDave et al. 
20101) . In this scenario, SMGs are massive galaxies sitting at the centers of deep poten- 
tial wells and fed by smooth accretion. 

GN20 is an interesting SMG in this regard. Situated at z=4.0 5, CO observations 
hint at a large spatial extent and ordered rotation (ICarilli et al.ll2010h . However, the pre- 
vious VLA observations had some severe spectral limitations; in particular, the limited 
bandwidth truncated the line profiles, and the reliance on continuum mode precluded 
any information on kinematics. We therefore obtained 124 hours on the Expanded Very 
Large Array (EVLA) to image the CO(2-l) emission in the GN20 field in the B- and 
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Figure 1 . Left: CO(2-l) velocity-averaged map over 780 km s _1 at 0. 19" resolu- 
tion. The circles identify molecular gas clumps and are the same size as the beam. 
The rms noise is 19 fi]y bm -1 , and contours are given in levels of Icr starting at ±2<x. 
Right: Example spectrum for Clump #2 at 39 km s _1 resolution. Velocity offset is 
with respect to the systemic velocity of GN20. 



D-configurations. We centered the two 128 MHz IFs at 45.592 GHz and 45.720 GHz, 
for a total bandwidth of 246 MHz or 1600 km s -1 . Each of the two IFs had 64 channels, 
corresponding to ~13 km s _1 per channel. The ~20 hours o f lower-reso l ution D -array 
data were presented in the EVLA Special Issue of ApJ by ICarilli et"aD (|201ll ). Here, 
we present the full dataset (B+D-configurations) on GN20, providing greatly improved 
spatial resolution and image fidelity. 



2. RESULTS & ANALYSIS 
2.1. Molecular Gas Clumps 

A velocity-averaged map is shown in Figured] (left). The gas appears to be resolved 
into multiple star-forming clumps, where we define as a clump anything with an H2 
column density greater than 6290 x (a/0.8) M Q pc~ 2 , or ~35 times the a verage surface 
brightness of local giant molecular clouds (GMCs: [Solomon et al.lll987l) . Taking the 
clump size as the beamsize, this corresponds to a velocity-averaged CO flux density 
threshold of 76 juJy beam" 1 , or the 4cr contour in the map. The clumps defined in this 
way have assumed sizes of 1.3 kpc (our beamsize), many orders of magnitude larger 
than even the giant molecular clouds in the nearby universe. 

Using this definition of a molecular gas clump, we identified seven clumps in 
GN20 (identified by the circles in Figure [TJ). We then extracted a spectrum at the peak 
position of each clump. While many of the clump spectra have low S/N, they allow us 
to estimate the properties of individual gas clumps for the first time. 

Using Gaussian fits to the spectra, we find that the clumps have linewidths of 
-100-500 km s" 1 (FWHM) and mass surface densities of 2,700-6,200 x (a/0.8) M Q 
pc~ 2 . Derived brightness temperature s range from ~11 K to 23 K (rest frame), ap- 



proaching the dust temperature (33K: iMagdis et aLll201 ll) and indicating that we are 



close to resolving the clumps. We find H2 masses of 3.6-8.2 x 10 9 x (a/0.8) M©. We 
therefore find that each clump makes up a few percent of the total mass of molecular 
gas. The combined mass of the clumps amounts to ~30% of the total gas mass. 
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Figure 2. Left: CO(2-l) 0th moment (i.e. integrated intensity) map at 0.19" res- 
olution. Contours start at (and are in steps of) 594 /vJy km s . Right: 0th moment 
map overlaid on the HST+ACS z850-band image. 



2.2. Dynamical Analysis 

For the dynamical modeling, we used the GALMOD task in the GIPSY package. We 
carefully masked the data to create 0th and 1st moment maps for input to the modeling. 
In the 0th moment map (Figure O, we recover more diffuse emission than was seen in 
the original velocity-averaged map of Figure [TJ We estimate a radius for the disk of 
1" ± 0.3", equivalent to ~7 kpc (± 2 kpc) at z = 4.055. The more stringent masking 
applied to the velocity field is evident in the 1st moment map (Figure [3l left), which 
shows a zoomed-in region of the 0th moment map. Although some noise is still present 
in the outskirts, a clear velocity gradient is apparent across the disk. 

By comparing different models to the data, we find that the best-fit model is a 
rotating disk with an inclination of i = 30° + 15°, a maximum rotational velocity of 
Vmax - 575 ± 100 km s , and a dispersion of 8 - 100 ± 30 km s _1 (Figure [3j right). 
While the resolution of the data make it difficult to constrain the exact shape of the 
rotation curve, we find that the velocity field is fully consistent with a rotating disk with 
a steeply-rising rotation curve that quickly flattens. We derive a dynamical mass for 
GN20 of 5.4 ± 2.4 x 10 11 M©. We will use this estimate to set limits on the CO-to-H 2 
mass conversion factor in the following. 

If we assume that GN20 is composed entirely of molecular gas, then we derive 
a conversion factor of 3.4 M© (K km s" 1 pc 2 )" 1 . If we instead assume roughly equal 
parts gas and stars, we derive a < 1.7 M© (K km s" 1 pc 2 )^ 1 . Here we have ignored 
the contributions from dust and dark matter, which will lower the limit even further. 
While a galactic conversion factor of ~4.3 M© (K km s" 1 pc 2 ) -1 may therefore still be 
possible (within the uncertainties), we consider this an extreme case. Our constraints on 
a are co nsistent with the lim it a < 1.0 derived using the local M gas IM c i ust vs. metallicity 
relation (Ma edis etal1l201il) . 



3. DISCUSSION 

The high-quality, high-resolution data on GN20 presented in this work reveal an ex- 
tended gas reservoir, 14 ± 4 kpc in diameter, which is consistent with a rotating disk. 
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Figure 3. Left: CO(2-l) 1st moment map at 0.19" resolution. Contours show 
steps of 100 km s _1 , with the contour near the green representing the systemic ve- 
locity. Right: Major axis position-velocity diagram for CO(2-l), taken at a position 
angle of 25°. The velocities on the vertical axis are relative. Greyscale and thin 
countours show the observed data, and thick contours show the best fit model. 

The data further resolve the reservoir into multiple, kpc-sized gas clumps which each 
constitute a few percent of the total gas mass and together comprise only 30% of the 
total gas mass. These observations may favor a scenario where the star formation is 
fueled by a process other than a major merger, e.g., cold mode accretion. To determine 
if this is in fact the case, a crucial question must be addressed: is it possible to have a 
highly-obscured starburst without a major merger? This is a question that, due to the 
faintness of this particular source, may have to await the launch of JWST 
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